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We studied the cosmological constraints on the Galileon gravity obtained from observa- 
tional data of the growth rate of matter density perturbations, the supernovae la (SN la), the 
cosmic microwave background (CMB), and baryon acoustic oscillations (BAO). For the same 
value of the energy density parameter of matter f2m,o, the growth rate / in Galileon models 
is enhanced, relative to the ylCDM case, because of an increase in Newton's constant. The 
smaller fim,o is, the more growth rate is suppressed. Therefore, the best fit value of fimfi 
in the Galileon model, based only the growth rate data, is quite small. This is incompatible 
with the value of fim,o obtained from the combination of SN la, CMB, and BAO data. On 
the other hand, in the ACDM model, the values of fimfi obtained from difi^erent observa- 
tional data sets are consistent. In the analysis of this paper, we found that the Galileon 
model is less compatible with observations than the ylCDM model. This result seems to be 
qualitatively the same in most of the generalized Galileon models in which Newton's constant 
is enhanced. 

Subject Index: 401, 453 

§1. INTRODUCTION 

The modifications to gravity at cosmological distances have recently received 
much attention as a possible explanation for cosmic acceleration. Many modified- 
gravity approaches have been proposed, including f{R) gravity [1])], scalar-tensor 
theories [2])-|l}], Gauss-Bonnet gravity [5D,[S])]) the Dvali-Gabadazde-Porrati (DGP) 
braneworld model [7j)-[9])], and Galileon gravity [TO |) - [M|) ]. 

Modified-gravity models need to be constructed so that Newton gravity is re- 
covered at short distances to ensure consistency with solar system experiments; the 
deviation from General Relativity (GR) can be allowed at large distances. 

There are two ways to recover General Relativistic behavior in the high density 
regimes relevant to solar system experiments. The first is the so-called chameleon 
mechanism [25p ] in which the field mass is different, depending on the matter density 
in the surrounding environment. If the field is sufficiently heavy in the regions of 
high density, a spherically symmetrical body can have a thin-shell around its surface 
so that the effective coupling between the field and matter is suppressed outside 
the body. This mechanism can be at work in viable f{R) gravity and scalar-tensor 
theory. 

The second method of recovering GR in regions of high density is to introduce 
non-linear self-interactions in the form of a scalar field (j). The simplest term of 



*•* E-mail: k_hirano@tsuru.ac.jp 



typeset using PTPTeX.cIs (Ver.0.9) 



2 



K. Hirano, Z. Komiya and H. Shirai 



the self-interaction is {'V(p)'^\I\(j). There have been attempts to restrict the form 
of the Lagrangian by imposing "Gahlean" symmetry, d^4> — >• d^4> + 6^ HH])]- The 
self-interaction of the form (V0)^n(/) respects Galilean symmetry in the Minkowski 
background. The self-interaction term induces decoupling of the field (p from gravity 
at small scales by the so-called Vainshtein mechanism [26])]. This allows the theory 
to recover general relativity for length scales smaller than the so-called Vainshtein 
radius [26])], which is consistent with solar system experiments (see Ref. [27|) for the 
implementation of the Vainshtein mechanism in these models). This term appears 
in the 4-dimensional effective theory of the DGP model. Unfortunately, the DGP 
model is plagued by the ghost problem [28]) ] and is incompatible with cosmological 
observations [2U]) ]. 

The presence of a nonlinear self-interaction of the form ^(0)(V0)^n(/), where ^{(f)) 
is a function of a scalar field (p, as a part of Galileon gravity models has recently been 
studied [I0 p - I24p ]. If (^((p) is constant, then this term respects Galilean symmetry 
in the Minkowski space-time, which keeps the equation of motion a second-order 
differential equation. This prevents the theory from introducing a new degree of 
freedom; perturbation of the theory does not lead to ghost or instability problems. 
Depending on the form of the function, some models do not have this "Galilean" 
symmetry. However, it is still possible to retain the desired properties. The equation 
of motion for the scalar field can remain of second order, which prevents the theory 
from having additional degrees of freedom. 

Galileon gravity can induce self-accelerated expansion of the late-time Universe. 
Therefore, inflation models inspired by Galileon theory have been proposed [30 p -l34 p ]. 
The evolution of matter density perturbations and conditions for the avoidance of 
ghosts and instabilities in Galileon models have also been studied [1^416]) . [T8]) . [T9]) . 
[M ]).[35]) ]. In Ref. ]36]), it is demonstrated that the integrated Sachs- Wolfe effect gives 
a stringent constraint on a subclass of the galileon model. 

In this paper, we study the cosmological constraints using observational data 
of the growth rate of matter density perturbations, as well as supernovae la (SN 
la) [37])]. the cosmic microwave background (CMB) [38])]. and baryon acoustic oscil- 
lations (BAG) [39p ]. for the action ()2Tp given below. In low-energy effective string 
theory there is a scalar field (p called the dilaton (that is coupled to the Ricci scalar 
R) of form F{(p)R |40p ]. Furthermore, the higher-order string correction contains a 
non-canonical kinetic term of form (Vcp)^ as well as a field self-interaction of form 
^{(p){'V(p)^\I\(p in the action. We accommodate non-linear field derivative terms using 
a Lagrangian of form K{(p, X) - G{(p, X)n(p [22) ]. 

We compare observed data of the growth rate of matter density perturbations 
to theoretical predictions. Whereas the background expansion history in modified 
gravity is nearly identical to that of dark energy models, the evolution of matter 
density perturbations in modified gravity is different from that of dark energy models 
(the cosmological constant is the standard candidate for dark energy). Thus, it 
is important to study the growth history of perturbations to distinguish modified 
gravity from models based on the cosmological constant or dark energy. Therefore, 
we computed the growth rate of matter density perturbations in Galileon cosmology 
and compared it to observational data. In addition, we studied the cosmological 
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constraints on the Galileon model obtained from SN la, CMB anisotropies, and 
BAO observations so as to test the vahdity of the model. 

This paper is organized as follows. In the next section, we describe the back- 
ground evolution and the evolution of linear perturbations in Galileon cosmology. 
In Sec. [3] we study the cosmological constraints on Galileon gravity obtained from 
observational data of the growth rate of matter density perturbations, SN la, CMB, 
and BAO. Finally, conclusions are given in Sec. |H 

§2. MODEL 

2.1. Background evolution 

The action we consider is of the form 



S = I d'^xy/—g 



-^F{(P)R + K{^, X) - G{<P, X)a^ + Lr, 



(2-1) 



where g is the determinant of the space-time metric Mp\ = (SvrG)^^/^ is the 
reduced Planck mass, G is the gravitational constant, R is the Ricci scalar, and (j) is a 
scalar field with a kinetic term X = — (1/2)(V0)'^. The function -F(0) depends on cp 
only, whereas K((j), X) and G{(j), X) are functions of both (p and X. Lm is the matter 
Lagrangian. We have used the notation: (V(^)^ = g^^^V ^cjN ycj), □(/> = g^^V ^Vv4>- 

Variation with respect to the metric produces the Einstein equations; variation 
with respect to the Galileon field (j) yields the equation of motion. For Friedmann- 
Robertson-Walker spacetime, the Einstein equations give 

SM^yFH'^ =p^ + p^- SM^yHF -K + 2XK^x + QH(j)XG,x - 2XG,<^, (2-2) 

- M^^FiSH^ + 2H) =pr + 2Ml^HF + M^^F + K- 2XG,x4> - 2XG,^, (2-3) 
and the equation of motion for the Galileon field gives 

{K^x + 2XK^xx + QH(j)G,x + 6H(j)XG ^xx — 2XG^(f)X — 2G^(f))4> 

H^HK^x + <pK,^x + 9H^G,x + 3if0G,x + GHXG^^x - 6HG^^ - G,^^<j))<p 

-K^ - GM^yH^F^^ - 3Mp2 iJF^ = 0, (2-4) 

where an overdot represents differentiation with respect to cosmic time t and H = 
a /a is the Hubble expansion rate. Note that we have used the following partial 
derivative notation: Kx = dK/dX and K^xx = d^K/dX"^ (and similarly for other 
variables), pm and pr are the energy densities of matter and radiation, respectively, 
and Pr is the pressure of the radiation. 

In this paper, as we are interested in the basic parameters of Galileon theory, 
we consider the following functions 

F{<P) = -^(P, (2-5) 

i^(</.,X) = 2^X, (2-6) 
G{4>,X) = 2i{4>)X, (2-7) 
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where uj is the Brans-Dicke parameter and ^(0) is a function of (j). This model 
is the Brans-Dicke theory with the fohowing self-interaction term: ^(0)(V0)^n(/). 

In this case, the Priedmann equation (|2-2|) and (j2-3|) can be written in the fol- 
lowing forms, respectively, 



3H^ 



{Pm + Pr + 



3H^ - 2H 



{Pr+P<p), 



(2-8) 

(2-9) 



where the effective dark energy density is defined as 

2 r ■ N / • \ 3' 



P<t> = 20 




and the effective pressure of dark energy is 

2 

a2/ 



20 





(2-10) 



(2-11) 



For the numerical analysis in this paper, we adopt a specific model in which 

^2 



(2-12) 



MpV2. 



where Tc is the crossover scale [TB]) ]. This coincides with the model by Refs. [T^ . [TC|) . 
Evolution of matter density perturbations for this model has been computed in Refs. 
11,[I5D. 

At early times, to recover general relativity, we set the initial condition ~ 
This reduces the Einstein equations ()2-8p and ()2-9p to their usual forms: 
3i72 = (p^ + /Or)/M2j and -Si?^ - 2H = Pr/M^j. This is the cosmological version 
of the Vainshtein effect, the method by which general relativity is recovered below a 
certain scale. At present, to describe the cosmic acceleration, the value of Vc must 
be fine-tuned. 

In Fig. [H we plot the effective equation of state in Galileon gravity, which 
is described as tt^gg = —p^/{^Hp^) — 1 using the effective dark energy density p^ in 
Eq. (j2-lUp . We set l7m,o = 0.30 for the Gahleon models in Fig. [H We can see the 
asymptotic behavior — )• — 1 for the future. A phantom- like behavior is found at 
low z. The smaller absolute value of the Brans-Dicke parameter (a; < 0) is, the 
smaller is at low z. The model in our paper is the Brans-Dicke theory extended 
by adding the following self-interaction term: ^{(f)){W (f))^\Z\(f). Thus, oj of this model 
is not exactly the same as the original Brans-Dicke parameter. If |w| is infinity with 
fixing the value of Om,o^ crossover scale Tc also approaches infinity, and the terms 
including ^(0) in Eqs. (j2-10p and (|2-lip behave non-trivially. In the case of fixing 
the value of i?m,0) even if |a;| is infinity, this model does not approach ylCDM. 
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Fig. 1. Effective equation of state in Galileon gravity as a function of redshift z for various 
values of the Brans-Dicke parameter uj. The parameters are given by J?m,o ~ 0.30. 



2.2. Density perturbations 

The evolution equation for the cold dark matter overdensity 6 in linear theory 
is governed by 

^ + 2//(^ - 47rGeff/3(5 ~ 0, (2-13) 

where G^s represents the effective Newton's constant in Galileon gravity. For the 
model specified by Eqs. ()2-5p . ()2-6p . and ()2-7p . the effective Newton's constant is 
given by 



G, 



1 



eflf 



1 + 



J 



2\2 



where 



J = 3 + 2w + 



2^ + 8^: 



(2-14) 
(2-15) 



The effective Newton's constant Gefj is close to Newton's constant G at early times, 
but increases at later times. 

We set the initial conditions 5 ~ a, 5 ~ d at early times. As we are interested in 
the difference between the growth of density perturbations in the Galileon gravity 
and that in the ylCDM model, we assume that initial conditions of matter density 
perturbations are same as the conventional ylCDM case. Solving the evolution equa- 
tion numerically, we obtain the growth factor 6/a for Galileon gravity. The linear 
growth rate is written as 

/ = ^^ (2-16) 

dm a 
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Fig. 2. Growth rate / in Galileon gravity as a function of redshift z for various values of the 
Brans-Dicke parameter lj. The parameters are given by i?m,o = 0.30. 
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Fig. 3. Growth rate / in Galileon gravity as a function of redshift z for various values of today's 
energy density parameter of matter Qm,o- The parameters are given by a; = —10000. 
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The growth rate can be parameterized by the growth index 7, as defined by 



The matter energy density parameter is defined as i7m = Pm/'^M^^H'^ . 

We plot the growth rate / in Galileon gravity in Fig. [2] and Fig. [3l a; is 
the (constant) Brans-Dicke parameter and f^mfl is the energy density parameter of 
matter at the present day. We set f2rn,o = 0.30 for the Galileon model in Fig. [2] and 
u = -10000 for the Gahleon model in Fig. [3l f2m,o for the ylCDM model is 0.30 in 
these figures. 

For the same value of f^m,o, the growth rate / in a Galileon model is enhanced 
compared to the ylCDM case, because of the enhancement of Newton's constant. 
The smaller i7m,o is, the more suppressed the growth rate is. We describe the 
observational data of Fig. [2] and Fig. [3] in the next section. 



In this section, we study the cosmological constraints on Galileon gravity so as to 
test the validity of the model. We carried out a detailed investigation of the allowed 
parameter region using the following observational data. 

3.1. Observational data 
3.1.1. Growth rate 

In Table H] we list the observational data used for the growth factor / from galaxy 
redshift distortions and Lyman-a forests, as compiled in Refs. HT ]) - t43p . adding the 
data from The WiggleZ Dark Energy Survey |44p ] from original references, provided 
to the right of the table. 

Table I. Currently available data for linear growth rates fobs used in our analysis, z is redshift; a 
is the la uncertainty of the growth rate data. 



(2-17) 



§3. OBSERVATIONAL CONSTRAINTS 



z 



fobs 



a 



Rcf. 



2.125 - 2.72 



2.8 - 3.6 
3 - 3.8 



2.2 - 3 
2.4 - 3.2 



2.6 - 3.4 



0.15 
0.35 
0.42 
0.55 
0.59 
0.77 
0.78 
1.4 
3.0 



0.49 
0.7 
0.73 
0.75 
0.75 
0.91 
0.71 
0.9 
1.46 
0.74 
0.99 
1.13 
1.66 
1.43 
1.3 



0.1 
0.18 
0.09 
0.18 
0.09 
0.36 
0.14 
0.24 
0.29 
0.24 
1.16 
1.07 
1.35 
1.34 

1.5 
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EH 
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The statistical function for the growth rate is defined as 

^f = T. — — ' (3-1) 

where a{z) is the la uncertainty in the fobs{z) data. 
3.1.2. Type la supernovae 

The Union2 SN la data [37|) ] consists of low z SN la data observed at the 
F.L. Whipple observatory of the Harvard-Smithsonian center for astrophysics [52])], 
intermediate z data observed during the first season of the Sloan Digital Sky Survey 
(SDSS)-II supernova survey |53p ]. and high z data from the Union compilation |54p ]. 
The Union2 SN la data used the SALT2 light curve fitter because it performs better 
than either SALT or MLCS2k2 in terms of the scatter around the best-fit luminosity 
distance relation [37])]. 

The 557 Union2 SN la data set [37]) ] gives the distance modulus at redshift 

l^obsi,Zi). 

The theoretical distance modulus (for a given model) is defined by 

;u(2;) = 5 logio^^L + /^o, (3-2) 
where is the Hubble-free luminosity distance given by 

"' = (' + ='1^0)"^' 

and /Uo is 

/xo = 5 logio + 25 = 42.38 - 5 logio ^, (3-4) 

where h is the Hubble constant Hq in units of 100 km s~^ Mpc^^. 
For the SN la data we have 



557 I- , s / \ n 2 

IJ,{Zi) - flobsiZi) 



i=l 



CTfiiZi) 



(3-5) 



where o"^ is the total uncertainty in the distance modulus. 

As the nuisance parameter /Uq (Eq. (j3-4p ) is model- independent, we analytically 
marginalize it as follows: 

xiNia = a-— , (3-6) 



where 

557 r . - 

(^l{zi) 



1^ ^27^^ ' i-^-O 



=1 "M^ 

557 



^-Z. ^27T^ ' (3-8) 
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and 

557 ^ 

3.1.3. Cosmic microwave background 

The positions of CMB acoustic peaks are affected by the expansion history of the 
Universe from the decouphng epoch up to today. To quantify the shift of acoustic 
peaks we use data points la, R, of Ref. 138j) (WMAP7), where la and R are two 
CMB shift parameters and 2* is the redshift at decouphng. For the flat Universe, we 
have 

R = ^/J2;^Hoda{z,), (3-10) 

where = Jq* dz/H{z) is the comoving angular diameter distance. 

The multipole la is defined by 

la = 7r^, (3-11) 



where rs{z^:) = dz/[H{z)y'3{l + 3 J7f,^o/ (4^^7,0(1 + z))}] is the sound horizon at 
the decouphng. Hi^ q and f2^^Q are today's density parameters of baryons and photons, 
respectively. 

For the redshift z^, there is a fitting formula by Hu and Sugiyama [55])]: 

z^ = 1048[1 + 0.00124(f2fe,o/i^)"°-^^^][l + 9i{^m,ohy'], (3-12) 

where 

0.0783(f?,,o/^^)-°-^38 

l + 39.5(I?b,o/i')°-^^=^' ^ ' 

- l + 21.1(I?,,o/i2)i.8i- (3-14) 

In our analysis of the Galileon model, we adopt the values h = 0.71 and fi^^Q = 
0.02258/1-2 [38])], and the standard value = 2.469 x IQ-^h-'^. 

For a flat prior, the 7-year WMAP data measured best-fit values are [38])] 









VCMB = 


R 


H 









302.09 ± 0.76 

1.725 ±0.018 I . (3-15) 
1091.3 ±0.91 

The corresponding inverse covariance matrix is [38])] 

2.305 29.698 -1.333 
CcMB = I 29.698 6825.270 -113.180 | . (3-16) 
-1.333 -113.180 3.414 

/ /a - 302.09 \ 
^CMB = i?- 1.725 , (3-17) 
\ z^- 1091.3 / 

and find the contribution of CMB is 

XcMB = -'^cmrC'cmb^cmb- (3-18) 



Thus, we define 
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3.1.4. Baryon acoustic oscillation 

For BAO analysis, we apply the maximum likelihood method using the data 
points of Ref. Ml (SDSS7): 

/ rs{za) - 0,1905 -^0,0061 \ 



Dv{0.35) 



0.1097 ±0.0036 



where rs{zd) is the sound horizon at the baryon drag epoch z^- For z^, we use the 
fitting formula by Eisenstein and Hu [56)]: 



where 
and 



bi = 0.313(I?^,o/i^)"°-^^^[l + Om7{nm,ohy-^^^], (3-21) 



62 = 0.238(I?^,o/i^)°-^^^- (3-22) 
The dilation scale Dy at the redshift z is 

nl/3 



Dviz) 



daiz 



H{z) 



(3-23) 



where da{z) is the comoving angular diameter distance. 
Thus, we construct 

rs{zi) _ 0,1905 



\ Di.(0.35) 

and using the inverse covariance matrix |39p ]. 



Xbao = ( ^rtS : ri 1 , (3-24) 



0.1097 



1 30124 17227 . 

^BAO ~ ivoov Qfinvv • 1,0 zo; 



'BAO \^ 17227 86977 
The contribution of BAO to is 

Xbao = ^baqC'bao^bao- (3-26) 

To determine the best value for, and the allowed region of, the parameters, we 
use the maximum likelihood method. 

3.2. Numerical results 

We now present our main results for the observational constraints on the Galileon 
gravity specified by Eqs. ([ZSD, ([22]), ([23), and (f2T2D . 

In Fig. m we plot the probability distribution of the energy density parameter of 
matter i7rra,o in the Galileon model from observational data of (individually) growth 
rate, SN la, CMB, and BAO, and also the combination of all data, where the other 
parameter is marginalized. 
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Fig. 4. ID probability distribution of the energy density parameter of matter i7m,o for the Gahleon 
model from observational data of Growth Rate, SN la, CMB, BAO, respectively, and the com- 
bination of All data. The line of CMB almost overlaps with the line of All. 



Considering only the growth rate data, the best fit value for the Galileon model 
is i?m,o = 0.166. This is considerably smaller than value obtained from observations 
of SN la, CMB, and BAO. Using only the growth rate data, we obtained the following 
constraint: 

f2m,o = 0.166tom9- (68.3% C.L.) (growth rate data). (3-27) 

By contrast, using the combination of growth rate, SN la, CMB, and BAO data, 
the best fit value for the Galileon model is f^mfl = 0.317; we obtained the stringent 
constraint as follows: 

l2rn,o = 0.317t[}:|][}^. (68.3% C.L.) (all data). (3-28) 

In the Galileon model, the bounds produced by each set of observations are incon- 
sistent. 

In Fig. [21 we plot the probability distribution of J7m,o for the ylCDM model 
from observational data of (individually) growth rate, SN la, CMB, and BAO, and 
also the combination of all data. 

In the ylCDM model, the bounds from each observational data set are consistent. 

In Fig. El we plot the probability distribution of the Brans-Dicke parameter oj 
for the Galileon model from observational data of growth rate, the combination of 
SN la, CMB, BAO data, and the combination of all data, where other parameter is 
marginalized. 
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Fig. 5. ID probability distribution of J?m,o for the TlCDM model from observational data of growth 
rate, SN la, CMB, and BAO, individually, and the combination of all data. The line of CMB 
almost overlaps with the line of All. 
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Fig. 6. ID probability distribution of absolute value of the Brans-Dicke parameter \u}\ (note that 
cj < 0) for the Galileon model from observational data of growth rate, the combination of SN 
la, CMB, BAO data, and the combination of all data. 
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Using only the growth rate data, we obtained the following constraint: 

u < -1060. (68.3% C.L.) (growth rate data). (3-29) 

Furthermore, from the combination of growth rate, SN la, CMB, and BAO data (all 
data), we obtained the stringent constraint as follows: 

uj < -36900. (68.3% C.L.) (ah data). (3-30) 

In Fig. [71 we plot the probability contom's in the {f2m,o, |u;|)-plane in the Galileon 
model. The contours show the la (68.3%) and 2a (95.0%) confidence limits, from 
the growth rate data, and from the combination of SN la, CMB, and BAO data. 
(|a;| is the absolute value of the Brans-Dicke parameter; u; < 0.) 
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Fig. 7. Probability contours in the (f?m,o, |tij|)-plane for the Galileon model. The contours show 
the la (68.3%) and 2a (95.0%) confidence limits, from the growth rate data, and from the 
combination of SN la, CMB, and BAO data, respectively. 

In the Galileon model, the allowed parameter region obtained using only the 
growth rate data does not overlap with the allowed parameter region obtained from 
the combination of SN la, CMB, and BAO data at all. 

In Table [III we list the best fit parameters, the values, and the differences 
of the Akaike information criteria (AIC) ^7)] and the Bayesian information criteria 
(BIC) [58j) ] for both the Galileon model and the ylCDM model. Here, we use the 
combination of growth rate, SN la, CMB, and BAO data (all data). The value, 
AIC change, and BIC change for the Galileon model are larger than those of the 
ylCDM model. In the Galileon model, as the bounds from growth rate data is 
inconsistent with the bounds from other observational data, the value for all 
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data is large. In the analysis of this paper, we find that the Galileon model is less 
compatible with observations than is the ylCDM model. This result seems to be 
qualitatively the same in most of the generalized Galileon models in which Newton's 
constant is enhanced. 

Table II. Results of observational tests using growth rate, SN la, CMB, and BAO data (all data). 
We consider the spatially fiat Universe. 



Model 


Best fit parameters 




ZiAIC 


Z\BIC 


ylCDM 


Omfi = 0.268 


550.674 


0.00 


0.00 


Galileon 


X2^,o = 0.316 


582.013 


33.339 


37.697 




oj = -2.69 X 10" 









In Fig. [HI we plot the growth rate / using the best fit models of the Galileon 
theory and ylCDM using only the growth rate data, and using the combination of 
SN la, CMB, and BAO data. 
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Fig. 8. Growth rate / in best fit models of Galileon theory and /ICDM from only the growth rate 
data, and from the combination of SN la, CMB, and BAO data. 

For the Galileon model obtained from the combination of SN la, CMB, and BAO 
data (excluding growth rate data), the growth rate of matter density perturbations 
is enhanced. On the other hand, in the ylCDM case. The model using only growth 
rate data is consistent with the model using the combination of SN la, CMB, and 
BAO data. 




Galileon from SN la + CMB + BAO 
ACDM form Growth Rate 
ACDM from SN la + CMB + BAO 
Observations 
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§4. CONCLUSIONS 

For the same value of fimfi-, the growth rate / in Gahleon models is enhanced 
compared to that of the ylCDM case, because of the enhancement of Newton's con- 
stant. The smaller l7m,o is, the more suppressed growth rate is. Therefore, the best 
fit value of i7m,o in the Galileon model using only the growth rate data is Qmfl = 
0.166, a value that is relatively small. 

By contrast, the value of i7m,o in the Galileon model obtained from the combi- 
nation of SN la, CMB, and BAO data is larger than that in the ylCDM case. This 
result is consistent with Refs. [23]) In Refs. [23]).[2l|) . the cosmological constraints 
using observational data of SN la, CMB, and BAO have been studied, but there the 
models are not the same as the one in our analysis. In this paper, we add the impor- 
tant constraints obtained from the growth rate data of matter density perturbations, 
compare the allowed parameter region from the growth rate data to the region from 
the combination of SN la, CMB, and BAO data. 

As seen in Fig. [71 in the Galileon model, the bound based only on the growth 
rate data is completely incompatible with the bound based on the combination of 
SN la, CMB, and BAO data. This result seems to be qualitatively the same in most 
of the generalized Galileon models in which Newton's constant is enhanced. On the 
other hand, in the /ICDM model, the bounds based on each observation data set are 
consistent. 

We also find that the upper limit of the Brans-Dicke parameter uj for the Galileon 
model is not constrained tightly only by growth rate data, but that the stringent 
upper limit of a; is obtained from the combination of growth rate, SN la, CMB, and 
BAO data. 

Galileon gravity is a fascinating model that can induce the accelerated expansion 
of the present Universe, meets the requirements of solar system experiments, and can 
avoid ghost and instability problems. However, through the analysis presented in 
this paper, we found that the Galileon model is less compatible with observations 
than the ylCDM model is. 

As the constraint from (only) growth rate data is not tight, more and bet- 
ter growth rate data are required to distinguish between dark energy and modified 
gravity. 
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